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Abstract 

In order to settle the fundamental question whether the nuclear forces involve 
the long range components, the S-wave amplitude of the proton-proton scattering is 
analysed in search for the extra singularity at z/ = 0, which corresponds to the long 
range force. To facilitate the search, a function, which is free from the singularities in 
the neighborhood of z/ = when all the interactions are short range, is constructed. 
The calculation of such a function from the phase shift data reveals a sharp cusp 
at 1/ = in contradiction to the meson theory of the nuclear force. The type of 
the extra singularity at z/ = is close to what is expected in the case of the strong 
van der Waals interaction. Physical meanings of the long range force in the nuclear 
force are discussed. Low energy p-p experiments to confirm directly the strong long 
range interaction are also proposed, in which the characteristic interference patten 
of the Coulonb and the Van der Waals forces is predicted. 
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1 Introduction 



In the Yukawa model of the nuclear force[|T|, the pi-meson plays the central roll, 
since the nuclear potentials are assumed to arise from the exchanges of a pion and 
a set of pions at least outside of the inner core region. Although the two-pion 
exchange potentials constructed directly from the meson theory are not successful 
in understanding the details of the data of the nucleon-nucleon scatterings, various 
phenomenological nuclear potentials have been proposed to reproduce the existing 
data. In the constructions of the nuclear potentials, two features are commonly 
shared by all the nuclear potentials, which are inherited from the meson theory of 
the nuclear force. One is the one-pion exchange (OPE) part of the potential and 
the other is that the range of the non-OPE part is shorter than l/(2/i). In terms of 
the singularities of the scattering amplitudes A{s, t), such features correspond to the 
one-pion exchange pole at t = /x^ and the cut of the two-pion exchange which starts 
at t = These analytic structures are common to all the amplitudes calculated 
from the various proposed nuclear potentials. 

On the other hand, in the nineteen sixties our views on the nucleons changed from 
'elementary particles' to composite particles, and moreover in the most important 
models of hadron, such as the QCD and the dyon model[0], the constructive forces 
of the composite particles are the Coulombic types. Therefore we cannot exclude 
the possibility for the induced long range force to appear between the composite 
particles. In particular, in the dyon model of hadron, the appearance of the strong 
van der Waals interaction between hadrons is a natural cosequence, because in such 
a model the hadrons are regarded as the 'magnetic atoms'!^. In general the long 
range force gives rise to a singularity at t = in the scattering amplitude A{s, t), and 
therefore the analytic structure is completely different from the case of the purely 
short range interactions. 

Since there is no apriori reason to believe that the strong interaction is a synonym 
for the short range interaction, it is desirable to search for the possible long range 
components of the nuclear force, whenever new precise data in \t\ <C 4/^^ come 
to be available. It is the difference of the analytic structure that allows us to do 
the search in the clear-cut way without being disturbed by the ambiguities such as 
different parametrization of the nuclear potentials in the inner region. Therefore, 
our aim in the present article is to search for the possible extra singularity of A{s, t) 
at t = 0, which is characteristic of the long range force, or to find the corresponding 
extra singularity in the partial wave amplitude. It is important that the search of 
the extra singularity can be carried out independently of the uncertainties of the 
spectra related to other singularities. The reason of such possibility comes from the 
fact that the location of the extra singularity is the end point of the physical region 
— 4z/ < t < 0, where the experimental data are available. Therefore, in contrast to 
the search for the short range force, our search of the extra singularity at t = does 
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not require any analytic continuation from the physical region. In our approach 
based on the analytic structure, the dispersion technique will turn out to be useful. 

In the search of the long range force, the required accuracy of the input data 
depends heavily on the power 7 of the threshold behavior of the spectral function 
At{s,t) at t = 0, where for small u and \t\ the power 7 is introduced by At{s,t) = 
ttC'P + ■ ■ ■. Since 7 and another power a, which appears in the asymptotic behavior 
of the long range potential as V{r) ~ —C/r°', are related by a = 27 + 3, our 
observation of the extra singularity in the amplitude becomes more and more difficult 
and requires the higher precision of the input data, as a increases. For example, the 
van der Waals potentials of the London type (a = 6) and the Casimir-Polder type 
{a = 7) imply 7 = 1.5 and 7 = 2.0 respectively, and therefore to observe these van 
der Waals forces is to recognize the sing ularities of the type C"(-t)3/2 (j,^2 log(-t) 
on the smooth back ground function of A{s,t). In order to observe the long range 
force, it is essential to obtain the high precision data at least in the region of small 

when the power a of the asymptotic potential is not small. 

In the hadron physics, the accuracies of the phase shift data of the S-wave proton- 
proton scattering in the low energy region are exceptinally high[0, and so it is 
better to analyse the S-wave amplitude /io('^) instead of A{s,t) in our search for 
the long range force. After the partial wave projection, the singularity of A{s, t) 
at t = ti becomes the left hand singularity at z/ = — ti/4 in hf){i>). Therefore the 
extra singularity at t = appears at z/ = 0, whereas the OPE pole changes to the 
logarithmic cut starting at z/ = — and the two-pion exchange spectrum starts 
at z/ = — /i^ in holu). In addition to these singularities, /lo(^) has the unitarity cut 
in z/ > 0. Moreover since we are considering the proton-proton scattering, hQ{i>) 
has also the cuts of the Coulombic interaction in — 00 < z/ < and of the vacuum 
polarization in —00 < z/ < — . 

In section 2, we shall construct a function KqIu) whose non-OPE part is free from 
the singularities in |z/| < /i^, when the forces are the short range types of the meson 
theory plus the electromagnetic interaction. If we take advantage that the scattering 
length is known within 0.05%, we can even consider the once subtracted function 
^once^^^ _ (^Koi^u) — ii'o(0))/z/. The merit to use the once subtracted function 
instead of KqIu) is that, since the power of the threshold behavior at z/ = changes 
from 7 to (7 — 1), the extra singularity becomes much easier to be obsreved, when it 
exists. In section 3, by using the phase shift data we shall evaluate numerically the 
once subtracted Kantor amplitude Kq'^'^'^{i>), and its non-OPE part Kq"''^'^{i') will be 
tabulated. The result of the evaluation is that Kq^'^'^{v) has a sharp cusp at z/ = 0. 

In sections 4 and 5, Kq^'^^{v) is fitted by using the spectrum of the long range 
force and that of the short range force respectively. It turns out in section 4 that 
the chi-square minimum occurs at 7 = 1.48 ~ 1.62, which is close to the 7 of the 
van der Waals interaction of the London type. In section 5, it is shown that the 
conventional spectrum of the short range interaction cannot reproduce the cusp. 
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In section 6, by assuming that the well-known mechanism to produce the van 
der Waals force is working, the lower bound of the strength of the underlying 
Coulombic force is estimated by using the inequality of the strength Ce of the Van 
der Waals potential. The lower bound of becomes 3.3 or 14 depending on the 
value of the radius of the composite particle, which comes from the measurement of 
the nucleon form factor or the distance where the deviation from the asymptotic form 
of the van der Waals potential becomes appreciable, respectively. In section 6, it is 
also pointed out that the Coulombic interaction between the magnetic monopoles is 
an imprtant candidate of the underlying super-strong Coulombic force, and the dyon 
model of hadron is explained briefly. Section 7 is used for remarks and comments, 
in which low energy experiments to confirm the long range force, by observing the 
characteristic destructive interference pattern of the Coulomb and the strong Van 
der Waals forces, are proposed . 

2 Selection of a Regular Function 

The low energy proton-proton scattering is prominent in its accuracy of the mea- 
surements in the hadron physics. Especially the S-wave amplitude /lo(i^) in the 
low energy region provides the ideal place to answer to the fundamental question 
whether the hadron-hadron interactions are short range except for the electromag- 
netic components. In the following investigations, we shall make use of the difference 
of the analytic structures of the amplitudes ho{i'). In general, for the short range 
interaction which arises from the exchange of a state with mass ^/tl, a singularity 
appears at u — —ti/A in the partial wave amplitude ho{i'). On the other hand, the 
long range interaction gives rise to a singularity at z/ = 0. Since the location of the 
extra singularity due to the long range force is the end point of the physical region 
1/ > 0, where the experimental data are available, we can expect to observe the 
singularity directly, if it exists, without making any analytic continuation from the 
physical region. The first thing we have to do is to construct a function which is 
regular at z/ = 0, when all the interactions are short range. Next thing is to eliminate 
the known near-by singularities from the function. The function with such a wider 
domain of analyticity serves to expose the extra singularity at z/ = 0, and makes it 
easier for us to observe the long range interaction, when it exists. 

Although the main aim of this section is to construct such an analytic function for 
the proton-proton scattering where the effects of the vaccum polarization as well as 
the Coulomb interaction are not negligible, it is instructive to start by constructing 
such a function for the neutron-neutron scattering first. It is well-known that the 
S-wave amplitude /lo(i^) has the unitarity cut in u > with the spectral function 
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Imhoi^u), where /lo(z^) relates to the phase shift by 



The most famous function, which is analytic at u = and therefore accepts the 
Taylor expansion of u when all forces are short range, is the effective range function 
Xo(z/) of Bethe. As it is well-known, Xq^u) is defined by 

Xo{iy) = V^cot6o{iy). (2) 

From Eqs.(l) and (2), the relation between /lo(^) and Xo(z/) is 



Since the one-pion exchange (OPE) contribution is 

the values of the coupling constant g'^/'in and the neutral pion mass fi are sufficient 
to eliminate the OPE cut from holu). In fact holu) = holu) — hl^iy) does not have 
the OPE cut. However Eq.(3) indicates that in order to eliminate the OPE cut from 
Xq^v) information on HehQ^u) in z/ < — as well as g'^/'iir and /i are necessary. 
Therefore the effective range function Xo(z/) is not adequate for our purpose to 
construct a function with wider domain of analyticity. 

Another possibility is the Kantor amplitude]^, which is defined by 

Ko{u) = ho{u) / du — . (5) 

TT JO V — V 

The Kantor amplitude does not have the unitarity cut, and the values of Kq{v) 
for real positive v can be evaluated directly from the experimental data, although 
the integration must be regarded as the principal value integration of Cauchy. It is 
straightforward to eliminate the OPE cut from the Kantor amplitude, and which is 
achieved by introducing the non-OPE Kantor amplitude Kq{v) by 

k,{u) = K,{u)-Kl''{u) , (6) 

where the OPE part of the Kantor amplitude is common to that of the partial wave 
amplitude, namely Kl'^{u) = h]^{v). It is useful to rewrite the Kantor amplitude 
introduced in Eq.(5) into the form of the contour integration. 
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where the closed contour C is shown in Figure 1. The merit to write Kq{v) in the 
form of Eq.(7) is that because of Eq.(3), when the effective range function Xq{v) is 
a polynomial of v or more generally a meromorphic function of y^, then the only 
singularities of hQ{v) on the first sheet of v are poles.Q Therefore by shrinking the 
contour C to a point, Kq{v) becomes the sum of the contributions of the poles on 
the first sheet of v or in the upper half AjA'-plane. In this way, we need not carry 
out the principal value integration of the rapidly changing function of Eq.(5). This 
technique to change the integration into the sum of the contributions of poles will 
be used in the actual calculation of Kq{v) in the next section. We are now in the 
position to examine Kq^u) in search for the extra singularity at = 0, if the precise 
data of the neutron-neutron scatterings were available. 
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Figure 1: The contour C of the integration in Eq.(7). vi and V2 are real poles on negative real 
axis, whereas Vc and v* are a pair of complex poles on the j/'-planc. Their values are listed in 
Table 2. 

Let us turn to the proton-proton scattering, where the vacuum polarization as 
well as the Coulombic interactions are important. The Kantor amplitude introduced 
in Eq.(5) does not satisfy our requirement of the analyticity. This is because the 
Coulombic cut in — oo < z/ < and the cut of the vaccum polarization in —oo < 
V < —ml still remain in the Kantor amplitude Ko{i>) of Eq.(5). The difficulties 
are by-passed if we remember the modified effective range function Xo(z/) of the 
proton-proton scattering, which is regular at u = and accepts the effective range 

^Strictly speaking V + v of the numerator of Eq.(3) causes a branch point at i/ = ^m? . 
Since v — ~m? is the far-away singularity, it does not affect our search for the extra singularity 
at 1/ = 0. However if we want to be more precise, we may redefine the amplitude by multiplying 
m/\/m? + V to hQ(v), and make necessary changes in the following calculations. 
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expansion, when all the forces are short range except for the terms of the Coulomb 
and of the vacuum polarization. The modified effective range function Xo(z/) for the 
phase shift S^i^u) is 

Xoiu) = ^^{il + Xo)cotS^-ta.nTo} + me^hiri)+me^ioir]) ■ (8) 

1-00 

In Eq.(8), two well-known functions with the Coulombic order of magnitudes appear, 
they are expressed using a new variable t] = me^/(2y^) : 

27111 °° 1 

In Eq.(7) tq is the phase shift due to the vacuum polarization potentialp 
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where irie is the mass of the electron and A = 2e^/37r = 1.549 x 10~^. Functions tq, 
Xoi 00 and ^0(^7) have the order of magnitudes of the vaccum polarization, and it is 
sufficient for our purpose to retain only the first term in A. In such approximation, 
they are expressed by -Fo('") and Go{r), which are the regular and the irregular 
Coulombic wave functions respectively, as 

r„ = -2„Ar*MW!:) (U) 

JO r 





£o(r]) = -A r d/-^[iCGo{r)r - iCGo{r))l^o] , (13) 
Jo r 

where /(r) is defined in Eq.(lO). Exact definitions of these functions are found in 
the paper by Heller [|]. 

By using the modified effective range function Xo(i^), we define the S-wave am- 
plitude ho^u) of the p-p scattering by 



The relation between /lo(i^) and the phase shift is obtained if we substitute X^iv) 
of Eq.(8) into Eq.(14), and which reduces to the well-known form 

h,{v) = ^^^(^p^e<^^U.^b^,{v) , (15) 
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if the functions related to the vacuum polarization are neglected. The form of ho{iy) 
of Eq.(15) is the same as that of the neutron-neutron scattering of Eq.(l) except for 
the factor Cq given in Eq.(9), which is the penetration factor. If we compare the S- 
wave amplitude of the p-p scattering with that of the n-n scattering, which are given 
in Eq.(14) and Eq.(3) respectively, a combination of functions {—me^h^q) — isJvC'^ 
appears in place of —i\fv. In order to investigate the analytic structure of holv), 
it is convenient to rewrite the combination as 

-me^^rj) -iy/uC^ ^ -iy/u + me^{\og{iri) -'ilj{l + irj)} . (16) 

Since the digamma function ip{z) has poles at non-positive integers, the poles on the 
ry-plane appear on the positive imaginary axis. In terms of which is me^/ {2r]) , 
the series of poles appear on the negative imaginary axisis and converge to ^/u = 0. 
It is the smallness of the fine structure constant and therefore of the residues 
of such poles that zeros of the denominator of Eq.(14) occur at points very close 
to the locations of the poles of {—me'^h{r]) — ii/z/Cg). Therefore the partial wave 
amplitude ho{i') of the p-p scattering has a series of poles on the second sheet of u, 
namely on the lower half plane of -y/i^, whereas on the first sheet of u the analytic 
structure of hoi^u) does not change compared to the case of the n-n scattering. This 
fact implies that the same definition of the Kantor amplitude Kq^u) introduced for 
the neutron-neutron scattering, which is given in Eq.(5), is valid also for the proton- 
proton scattering, as long as we evaluate lmho{u') of Eq.(5) from Eqs.(14) and (16). 
Therefore the Kantor amplitude of the p-p scattering Ko{i') constructed in this way 
is free from the singularities in the neighborhood of z/ = 0, and so does not have the 
cut of the vacuum polarization as well as that of the Coulomb interaction. 



3 Numerical Calculation of the Kantor 
Amplitude 

Since the scattering length of the S-wave amplitude of the p-p scattering is known 
with high precision, we shall analyse the once subtracted Kantor amplitude, which 
is 

Kr» - ^"^"^^'^^^^^ (17) 

Merits to use the once subtracted Kantor amplitude rather than Ko[i') are twofold. 
The first one is the higher convergence of the integration 

U TT Jo — u) 

namely the uncertainty arising from the lack of the very accurate data in the higher 
energy region is largely suppressed because of the extra factor v' in the denominator 
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of the integrand. The second one is the change of the power of the threshold behavior 
of the spectral function from C'{—v')"' to — C"(— z/')'''"^, and which makes it easier 
for us to recognize the extra singularity at = in Kq"'^^{v). Since the power a, 
which appears in the asymptotic behavior of the potential as V{r) ~ —Ca/r"'-, and 
the power 7 of the spectral function Im ho{h'') are related by 

a = 27 + 3 , (19) 

7 is equal to 3/2 for the van der Waals potential of the London type, and which is 
expected to occur in the dyon model of hadronp|0. Therefore in such a case the 
behavior of the once subtracted Kantor amplitude in the neighborhood of z/ = is 
that 

= const. + C'v^ + ■ ■ ■ with C >0 , (20) 

and we must observe the extra singularity as a cusp aX u = unless the coefficient 
C is very small. 

Let us represent the effective range function of the p-p scattering Xo(z^) by a 
meromorphic function of u: 

Xq{v)-co-civ V bo + biu 

1 ) = ao + aiu + -—— ■ — - , 21 

u Up do + di^ + ^ 

in which up is the location of the pole of Xo(z^). — Cq and Ci are the inverse of the 
scattering length and the effective range devided by 2 respectively. The r.h.s. of 
Eq.(21) involves six free parameters to be fitted. In the fitting, the Compton wave 
length of the neutral pion will be used as the unit of the length. In the search of 
the free parameters, cq, ci and up are fixed beforehand which are0 

Co = 0.18698 , ci = 0.9541 and vp = 6.39 . (22) 

Accuracies of Cq and Ci are as high as 0.05% and 0.6% respectively. On the other 
hand, uncertainty of z/p is around 3 per cent. However reasonable accuracy of up is 
sufficient for our purpose to search for the extra singularity at z/ = . 

The S-wave phase shifts of the p-p scattering by the Nijmegen group 
are used as the input to evaluate the effective range function Xq^u). The once 
subtracted form of Xo(z/), in which the pole at z/ = z/p is eliminated by multiplying 
a factor (1 — u/up), is displayed in figure 2. The points in the figure come from the 
multi-energy phase shifts of Nijmegen-90 p[. 

The box at the left lower corner in the figure is enlarged and displayed in figure 
3, and in which the low energy data in the energy range 0.1 MeV. < Tiab < 30 MeV. 
of Nijmegen-88 are shown with their error bars. The data in 0.5 MeV. < Tiab < 
100 MeV. are fitted by the meromorphic function of Eq.(21), and the six parameters 
are determined by the chi square search. The searched curve is shown in figures 2 
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Figure 2: (Xo(i^)-co-cit/)(l-i//j/p)/i/ of 
Eq.(21) is plotted against Tiat, where Xoi^) 
is the effective range function of the S-wave 
p-p scattering. The points are the data of 
Nijmegen 90. The box of the left lower corner 
is enlarged in fig. 3. 



0.01 



-0.01 



-0.02 



-0.03 




Figure 3: Enlarged graph of {Xo{i^) - cq - 
civ)[l — vlvp)lv with error bar is plotted 
against Tiah in Tiah < 30 MeK, in which the 
points are the data of Nijmegen 88 and 90. 
The parameters of the fitted curves are listed 
in Table 1. 



and 3 as the full line, and whose parameters are tabulated in the column [sc] in 
Table 1. On the other hand the dash curve in the figures is a different fit by a 
meromorphic function, in which the six parameters are fixed by the six data points 
at Tiah = 1, 5, 10, 25, 50, 100 MeV. of the Nijmegen-90|1. The parameters of the 
fixed curve are tabulated in the column [fc] in Table 1. 



Table 1: Parameters of the meromorphic form of the 
effective range function introduced in Eq.(21) 





[uc] 


[sc] 


[Ic] 


[fc] 




-0.2579 


-0.2703 


-0.2815 


-0.3295 


ai 


0.15278 


0.15504 


0.15694 


0.16328 


bo 


0.01009 


0.01983 


0.03463 


0.09272 


h 


0.1862 


0.2098 


0.2337 


0.39475 


do 


0.0379 


0.0720 


0.1218 


0.2782 


di 


0.812 


0.927 


1.061 


1.643 



In the same table, there are also columns of the curves of [Ic] and [uc], in which 
the parameters are determined by fitting to the lower fringe and to the upper fringe 
of the error bars in 0.5MeV. < Tiah < 100 MeV. respectively. Such parameters are 
necessary to estimate the errors of the Kantor amplitude. 
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It is remarkable that in figures 2 and 3, the curvature of the curve in the small 

T;„5 region increases very rapidly as Tiab decreases to zero. If we remember that the 
starting points of the spectra of the one-pion exchange(OPE) and of the two-pion 
exchange(TPE) are Tiab = —9.71 MeV. and —38.83 MeV. repectively, it is not easy 
to understand the rapid change of the curvature in the low energy region such as 
< Tiab < 10 MeV. . On the other hand, if we accept the long range force, the 
appearance of such a cusp at Tiab = is what is expected. Since the separation 
of the OPE term from Xo(z/) is not simple, in the present article we shall examine 
K^^'^^iy) because, contrary to X^iy), the separation of the OPE term from the 
Kantor amphtude is straightforward. 



Table 2: Poles and residues of hQ{u) for meromorphic Xo(z/) 





[uc] 


[sc] 


[Ic] 


[fc] 




-0.04951 


-0.08487 


-0.12894 


-0.18430 


n 


0.005345 


0.016313 


0.044157 


0.15756 




-0.54699 


-0.58917 


-0.63853 


-0.87572 


r2 


2.5629 


2.8719 


3.1755 


5.0604 


Re Uc 


-3.4550 


-2.8526 


-2.4303 


-1.3318 


Im Uc 


-6.0285 


-6.0191 


-5.9984 


-6.0552 


Re Tc 


-25.427 


-22.461 


-20.562 


-16.855 


Im Tc 


43.256 


36.108 


31.607 


21.008 



In Table 2, the locations of the poles of ho{i'), which is introduced in Eq.(14), 
on the first sheet of u and their residues are tabulated. The columns correspond 
to those in Table 1, and each of them has two real poles Ui and 1^2, and a pair of 
complex poles {uc, u*) respectively, on the first sheet of u. 

In Table 3, K^'^^{i'), which are the once subtracted Kantor amplitude of Eq.(17) 
minus the one-pion exchange contribution, are tabulated for [sc] and [fc]. Because 
of Eq.(7), 

Kr» = + + J + Ao - (OPEC) (23) 

Ui{u-Ui) U2{U-U2) Uc{u-Uc) 

, where OPEC is the one-pion exchange contribution. The correction term Aq arises 
from the approximation to Xo[i>) by the meromorphic function fitted to the phase 
shift in the energy range of 0.5 MeV. < Tiab < 100 MeV.. However the numerical 

value of Aq is extremely small, namely less than 2 x 10"'^ or less than 10% of 



10 



Table 3: Non-OPE part of the once subtracted Kantor amplitude Kq^^^{v) 





/77 

/ ly 


-Llab 






\ 


IV 


-Llab 








=1) 


(MeV.) 


\sc] 

L J 


\fc] 

L •/ J 


A T^once 
AAq 


(m 


= 1) 


(MeV.) 




\fc] 







05 


0.097 


-5.188 


-4.974 


2803 


1 


15 


51.35 


-2 


651 


-2.655 


0.0126 





10 


0.388 


-5.049 


-4.892 


2069 


1 


20 


55.92 


-2 


581 


-2.584 


0.0116 





15 


0.874 


-4.864 


-4.771 


0.1307 


1 


25 


60.67 


-2 


515 


-2.517 


0.0105 





20 


1.553 


-4.671 


-4.629 


0.0747 


1 


30 


65.62 


-2 


451 


-2.453 


0095 





25 


2.427 


-4.493 


-4.479 


0.0418 


1 


35 


70.77 


-2 


390 


-2.391 


0086 





30 


3.495 


-4.335 


-4.332 


0.0257 


1 


40 


76.11 


-2 


331 


-2.332 


0.0078 





35 


4.757 


-4.196 


-4.194 


0.0193 


1 


45 


81.64 


-2 


273 


-2.274 


0.0071 





40 


6.213 


-4.070 


-4.064 


0.0175 


1 


50 


87.37 


-2 


218 


-2.218 


0065 





45 


7.863 


-3.953 


-3.943 


0.0174 


1 


55 


93.29 


-2 


163 


-/.1d4 


0059 





50 


9.708 


-3.842 


-3.830 


0.0176 


1 


60 


99.40 


-2 


110 


-2.111 


0.0054 





55 


11.75 


-3.734 


-3.721 


0.0177 


1 


65 


105.7 


-2 


059 


-2.059 


0.0050 





60 


13.98 


-3.628 


-3.617 


0.0175 


1 


70 


112.2 


-2 


008 


-2.008 


0.0046 





65 


16.41 


-3.525 


-3.516 


0.0170 


1 


75 


118.9 


-1 


958 


-1.958 


0.0043 





70 


19.03 


-3.423 


-3.417 


0.0164 


1 


80 


125.8 


-1 


909 


-1.909 


0.0040 





75 


21.84 


-3.325 


-3.322 


0.0158 


1 


85 
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the error of K^'^^{i') for < z/ < 3. The smallness of Aq comes from the facts 
that firstly the meromorphic function reproduce Xo(z/) well in the wider domain of 
Tlab ~ ISOMeV., secondly in the neighborhood of Tiab = 248 MeV., where the phase 
shift passes through zero, the spectral function Im hQ{v')/v' is small and thirdly in 
higher v' region the spectral function ||10|| is suppressed by the factor z/', because we 



are computing the once sutracted Kantor amplitude. The last column of Table 3 is 
the error of K^"''^'^ (u) , which is defined by 

Infi gure 4, —Kq^'^'^^u) is plotted against Ti^b, and in which the diamonds are [sc] 
and the triangles are [fc] respectively. Figure 5 is the enlarged graph of the lower 
energy part of figure 4. The curves in the figures are that {L)^ is the 3-parameter 
fit by the spectral function of the long range force, whereas {sa)^ (dotted curve) 
and (56)3 (dash curve) are the fits by the short range potential with three free 
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parameters. Details of the curves will be explained in sections 5. A characteristic 



feature of —Ka 



is that it has a sharp cusp at z/ = 0, contrary to what is 



expected in the meson theory of the nuclear force, where the location of the two-pion 
exchange spectrum is 'far away' compared to the width of the cusp. In the one-pion 
exchange term of Eq.(6), (7^/47r 
However modifications to the different values | 



14.4 is used as the tt-N coupling constant [|n|][|l2|. 

I of g'^/4:TT are straightforward if 



we remember Eq.(4). It must be pointed out that as the coupling constant g /4:Tt 



decreases the cusp of —Kq^ 



becomes more prominent. 
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Figure 4: -K^'^^'^iv) is plotted against Tiab 
in Tiab < 125 MeV.. The diamonds and tri- 
angles are —Kq"'^'^{i') calculated from Xo{i') 
of the [sc] and the [fc] fits, whose parame- 
ters are listed in Table 1, respectively. The 
curve is the fit by the spectrum of the long 
range force in the range 0.6 MeV. < Tiab < 
125 MeV. . 

In the next two sections, Kn^'^'^iv] 



Figure 5: The enlarged graph of 
_Koncei^^-^ of [sc] is plotted in Tiab < 
20 MeV.. The three curves are 3- 
parameter fits to the data in 0.6 MeV. < 
Tiab < 125 MeV. The curve {L)^ is 
the fit by the spectrum of the long range 
force, whereas the curves of {s 0)3 (dots) 
and {s 6)3 (dash) are the fits by the short 
range potentials. 

of Table 3 will be analysed to understand 



the physical meanings of the cusp at u = 0. In section 4, the cusp will be fitted 
by the spectrum of the long range force, in which the power 7 of the threshold 
behavior of the spectrum will be searched. On the other hand, in section 5 it will 
be examined whether the short range nuclear potential V{r) can reproduce well the 
cusp of ^o°''''''('^)- 
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4 Fits by the Spectrum of the Long Range 
Interactions 

Figures 4 and 5 indicate that — A'q"^^(z/) has a cusp pointed upward aX u = 0. If there 
exists the spectrum of the long range force with the attractive sign, appearance of 
such a cusp in —Kq"-'^^{i') is a natural consequence. In this section, we shall determine 
the power 7 and the strength of the threshold behavior of the spectral function at 
z/ = 0. More conventional fit will be tried in the next section, in which the possibility 
to understand the cusp of — Kq'^'^'^ (u) in terms of the ordinary short range potential 
of the nuclear force will be considered. The integral representation of the scattering 
amplitude A{s,t) of the non-OPE part is 

where to and Uo are 4yU^ for the short range potential whereas zero for the long 
range interaction, respectively. The tilde such as A means the function in which 
the OPE part is already subtracted. Since we are considering the scattering of the 
identical fermi particles, the signs between the integrations of Eq. (25) are plus for 
the spin-singlet and minus for the spin-triplet states respectively. Moreover the 
spectral functions of t and u channels are the same, namely At{s, *) — A„(s, *). The 
partial wave projection of Eq.(25) gives 

1 r°° ~ 1 t 
ho{iy) = - dtA{s,t)—Qo{l + —) + {t^u). (26) 

TT Jo IV IV 

Therefore the spectral function of the left hand cut of the S-wave amplitude h^iy) 
is 

~ 1 /■— 41/ r—A,V 

lmho(u) = —{ dtAt{s,t)+ duAu{s,u)} for vK-t^/A. (27) 
41/ Jto Juo 

Let us search the value of the power 7 of the spectral function of At{4m'^,t) 
at t = or of ImhQ^u) at u = 0. Because of Eq.(27), the powers 7's of these 
functions are the same. Since ylj(4m^,t) directly relates to the potential by the 
Laplace transformation 

1 AT 

rV(r)^ dt'AMrn'^ty-'^ , (28) 

7rm^ Jo 

we shall assume the functinal form of 74j(4m^, t) and search the parameters involved. 
In the search of the power 7, three-parameter form of the spectral function will be 
used, which is 

it(4m^t) = 7rC"t^exp[-/5t] , (29) 
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where the exponential factor is necessary to make the integration Eq.(26) conver- 
gent .Q 

In the search, 68 points of K^^^{u) in O.GMeV. < Tiab < l25MeV. are fitted. 
The fitted points are chosen to be equispacing with respect to or more precisely 
= 0.025n with 5 < n < 72. The parameters (7 , /5 , C" ) and x of the best fits 

are 

7 = 1.543, /5 = 0.06264, C" = 0.1762 and x = 0.441 for [sc] (30) 
7 = 1.569, /5 = 0.06544, C" = 0.1734 and x = 0.343 for [fc] (31) 

respectively. The curves in figures 4 and 5 are the fits of [sc] (full line) and [fc] (dash 
line) respectively. 

In figure 6, xil) is plotted against 7, in which the full line is [sc] and the dash 
line is [fc] respectively, xil) is obtained by making the chi-square search for fixed 7. 
From the curves, minimum points are determined, which are 7mm = 1-543 ± 0.055 
and 1.569 ± 0.044 for [sc] and [fc] respectively. The range of the uncertainty of 7mm 
is estimated by the condition that xil) < ^■'^^Xmin, in which Xmin is the minimum 
value of x(7)- If is remarkable that the value of the power 7 is close to what is 
expected in the case of the van der Waals potential of the London type namely to 
7 = 1.50. 



gamma vs chi gamma vs chi 




gamma gamma 



Figure 6: xil) are plotted against 7. The 
curves are the fits to the —K^'^'^'^{v) of [sc] 
(full line) and of [fc] (dash line) respectively. 
The minimum is close to 7 = 1.5 of the Van 
der Waals force of the London type. 



Figure 7: xil) are plotted against 7 for 
[sc]. The three curves correspond to the 
different fitted regions. They are Tiab = 
(0.6 - 30)MeF. , (0.6 - l2h)MeV. and 
(0.6 ~ 175)MeV. respectively. 



However before we discuss about the physical meaning of the value that 7 ~ 1.5, 
it is necessary to examine the dependency of 7mm on the fitted range of energy. In 



^Moreover the parameter /3 plays the roll to incorporate the two-pion exchange spectrum in 
Eq.(29). 
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figure 7, three curves of xil) with different fitted energy domains are shown. They 
are (0.6 ~ 30) MeV. and (0.6 ~ 175) MeV., in addition to the domain (0.6 ~ 125) 
MeV., which was shown in fig. 6. The figure indicates that as the energy domain 
shrinks, 7mm moves slowly to samaller value. For 0.6Mel^. < Tiab < SOMeV. the 
best fit is that 7mm = 1.467 ±0.111 and Xmin = 0.423 for [sc]. On the other hand, as 
the energy region expands, 7mm increases and moreover the value of Xmin increases 
rapidly. Therefore the 3-parameter function of Eq.(29) is not suitable for the fit to 
so many data points. In fact, for the domain O.GMeV. < Tiab < 175MeV., the curve 
in the figure indicates that 7mm = 1.695 ± q^^q and Xmin = 1-410 for [sc]. Therefore 
even if we change the energy domain, the power 7 still remains in the neighborhood 
of 1.5 . 



5 Fits by the Spectral Function of the Short Range 
Interactions 

In the previous section, it was found that the cusp in KqIu) obtained in section 3 
was reproduced well by the spectral function At{^m?, t) of the long range interaction 
given in Eqs.(29),(30) and (31). However it is important to ask whether the cusp can 
be understood also in the framework of the more conventional approach, namely in 
terms of the spectral function of the short range interaction expected in the meson 
theory. It is well known that, in the meson theory of the nuclear force, the spectral 
function of the two-pion exchange is 

, in the neighborhood of the two-pion threshold at t = 4/i^. The power 1/2 of 
Eq.(32) comes from the functional form of the phase volume of the two-particle 
state exchanged. The coefficient 6*2^ can be estimated from the phase shifts of the 
TT-iV and the vr-vr scatterings by using the generalized unitarity of the dispersion 
techniqueP|. From Eqs.(27) and (32), the spectral function of the left hand cut of 
/lo(^), which is the same as that of i^'o(z/), becomes 

I^M^') = | ^C'l!,(-z.'-l)3/2 + ... In y<-\ . ^33) 

Since the power of the spectral function at the threshold u = —1 increases to 3/2, 
the distance of the average location of the spectrum from z/ = is much larger 
compared to that of the threshold point. With such a spectral function, it will be 
not easy to reproduce the cusp of Kq'^'^'^{i') observed in u <^ 1, where in our unit 
1/ = 1 corresponds to Tiab = 38.83MeV! . On the other hand, the updated nuclear 
potentials are believed to fit well to the precise phase shift data. [p!3| ll^ |1T5[| |[T6[] . 
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In this section, we shall examine the updated nuclear potentials more closely, 
and in which we shall concentrate on the cusp at z/ = 0. Among various potentials, 
we shall choose the regularized Reid potential updated by the Nijmegen group (Reid 
93) as an example ||15||, in which 50 phenomenological parameters are searched to fit 



to 1787 p-p data and 2514 n-p data in the energy range of < Tiab < 350MeV.. It 
is remarkable that by using such potential they achieved excellent fits = 1-03 per 
datum, which is close to the chi-square value 0.99 of the energy dependent phase 
shift analysis. In the Reid93, the S-wave potential of the non-OPE part of the p-p 
scattering involves five parameters: 

VppCSo) = A2Y{2) + A^YiS) + A,Y{4) + A,Y{5) + A,Y{6) , (34) 

where the regularized Yukawa functions of the range 1 / (p/x) are 

1 A 2 2 2 

Y{p) = p-[e-^-^^ - e-^^il + ^ Ar)] (35) 

and is the average masses of the three types of pions and A = 8/i. Contrary to the 
original Reid potentialfl^, this potential has a characteristic feature that the Y{2) 



term exists. Such term gives rise to a pole at t = Ajj? in the scattering amplitude 
A[s,t), whereas in the meson theory t = 4/i^ is merely the starting point of the 
continuous spectrum of the two-pion exchange. Therefore it is not desirable from 
the point of view of the meson theory to have the Y{2) term in Eq.(34). On the 
other hand the inclusion of the Y{2) term is necessary, in fitting the very precise 
low energy data of the p-p scattering, to reduce the the chi-square value. The 
contradictory situation already indicates the difficulty to reproduce the cusp from 
the far-away spectral function of the meson theory of the nuclear force. 

In figure 5, results of the 3-parameter fits to the data points in 0.6 MeV. < 
Tiab < 125 MeV. are compared. The curve {L)^ is the fit by the spectrum of the 
long range force, while the curves of (s 0)3 and (s 6)3 are the fits by the short range 
potentials of the Reid types of Eq.(34). The (s 0)3 is the sum of Y{2), ^(3) and 
y(4), on the other hand (s 6)3 is the sum of V(3), V(4) and 1^(5), where Y{p) is 
the regularized Yukawa potential of range 1/p defined in Eq.(35). The x- values per 
datum of (1^)3, (50)3 and (56)3 are 0.441, 1.82 and 3.11 respectively. It is evident 
that the short range potentials with three parameters cannot reproduce the cusp at 
V = {]. It is interesting that the inclusion the Y{2) term always reduce the value, 
this is because the lack of the spectrum of small t is the reason of the deviation from 
the data points of K^"-^^{iy). 

In figure 8, the seven curves are the fits to —Kq^'^'^{v) of [sc] in the off-cusp 
region Tiab = (20 ~ 125)MeV.. The fitted curves are then extrapolated to Tiab = 0. 
The (L)3 (full line) curve is the 3-parameter fit by the spectrum of the long range 
force of Eq.(29) as before, on the other hand (s a)„ (dotted line) and (s b)n (dash 
line) are the n-paremeter fits by the Reid type potentials. The differences of (s a)„ 
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5.2 




Figure 8: The seven curves are the fits to — _R'o"'^^('^) of [sc] i'^ t^^e ofF-cusp region Tiab = 
(20 ~ 125)MeV.. The fitted curves are then extrapolated to Tiat = 0. The {L)^ curve is the 
3-parameter fit by the spectrum of the long range force, on the other hand (s a)„ (dotted curves) 
and (s b)n (dash curves) are the n-paremeter fits by the Reid type potentials. The differences of 
(s a)n and (s 6)„ are that, although (s a)„ involves 1^(2) term, in (s 6)„ the Yukawa term of the 
longest range is Y{3). 

and (s b)n are that, although (s a)n involves Y{2) term, in (s 6)„ the Yukawa term 
of the longest range is Y{3). In general as the numbers of the parameters n increase, 
the curves come closer to the data points. However the curves of Fig. 8 indicate that 
the curves of the short range interaction [s b)n cannot reproduce the cusp of the 
data points at z/ = 0. 

In the meson theory, the nearest singularity of —Kq'^'^'^{i>) occurs at p = —1 or 
in terms of Tiab at —38.83 MeV., which is the threshold of the two-pion exchange 
spectrum. Since Eq. (33) indicates that the power of of the threshold behavior of the 
two-pion exchange spectrum at z/ = —1 is 3/2, the distance to the average location 
of the spectrum from the origin is much larger than one. What is important is 
that the length of our extrapolation (~ 0.5) is short compared not only to the size 
of the fitted domain (~ 3.0) but also to the distance to the left hand spectrum 
(> 1). Therefore if the analytic structure expected in the meson theory is correct, 
the extrapolated curve must approximately follow the data points. However figure 
8 indicates that this is not the case. On the other hand, the extrapolated curve of 
the long range force {L)^ traces the data point nicely. Therefore we conclude that 
the analytic structure of the amplitude of the nucleon-nucleon scattering is different 
from what is expected in the meson theory of the nuclear force, rather the nuclear 
force involves the large component of the long range force. 
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6 Strong van der Waals Interaction as a 
Candidate of the Nuclear Force 



In section 2, we introduced the S-wave Kantor amplitude of the proton-proton scat- 
tering KqIu), which was free from the singularities of the vacuum polarization and 
the Coulombic interaction as well as the unitarity cut. By subtracting the one-pion 
exchange contribution from Koiiy), the domain of analyticity expands further, and 
the nearest singularity of the non-OPE Kantor amplitude Kf){i>) occurs at z/ = — 1, 
and which is the branch point of the two-pion exchange spectrum. By using the 
phase shift data of the Nijmegen group^ [||, the once subtracted non-OPE Kantor 
amplitude K^'^^{i') defined in Eqs.(6)and (18) was calculated in section 3. The ap- 
pearance of the cusp at z/ = is rather surprising, because in the meson theory of 
the nuclear force, Kq^'^^{v) is regular at least in \v\ < 1. 

It turns out that the conventional fits of the spectra of the short range forces 
cannot reproduce the cusp properly even if the potentials with five parameters given 
in Eq.(34) are used in the fitting. On the other hand, if we use the spectra of the 
long range force, three parameters of Eq.(29) are sufficient to reproduce the cusp of 
Xonce^^y In summary, the short range condition, that Im/io(z^) = in z/ > —1, is too 
severe to reproduce the cusp in |z/| <^ 1. Therefore it is natural to regard the cusp as 
an effect of the long range interaction between hadrons. The power 7 of the threshold 
behavior was searched in section 4, and it turned out that 7 = 1.543±0.056 from the 
multi-energy data [sc] and 7 = 1.569±0.046 determined from the fit [fc], respectively. 
Since the power a of the asymptotic behavior of the long range potential V{r) is 
related to the power 7 by a = 27 + 3, the values of 7 imply a = 6.09 ±0.11 and 
a = 6.17 ± 0.09 for [sc] and [fc] respectively. 

It is remarkable that the values of a are close to that of the van der Waals 
potential of the London type ( a = 6 ). Although a can assume any value from 
— 1 (Coulomb) to infinity (Yukawa), Nature chooses a special value in the neigh- 
borhood of a = 6. Therefore it is difficult to believe the occurrence of the value 
a = 5.98 ~ 6.26 is merely accidental, rather the ordinary mechanism to produce 
the van der Waals interaction must be working. It is well-known in the atomic 
and molecular physics that, when the underlying Coulombic force exists, the dipole- 
dipole interaction occurs between the neutral composite particles, and the two-step 
process of such interactions causes the attractive potential [jTSf. The form of the 
potential between particle 1 and particle 2 is 

^ _^ (j^MEht , (36) 

n — ^0 

where H' is the hamiltonian of the dipole-dipole interaction 

H' = -^{^1^2 + ym - '2Z1Z2) , (37) 
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in which is the strength of the underlying Coulombic interaction. 

It is important to notice that each term in the summation of Eq.(36) is positive 
definite. Therefore if we replace all the denominators [En — Eq) in Eq.(36) by 
the first excitation energy (Aii^i), then by using the closure property the upper 
bound of the strength Cq of the van der Waals potential, which is introduced by 



V{R) = -Ce/R^ + ■ ■ -, is obtained. It isp 



< (AEi) - 3 (AEi) ■ ^^^^ 

On the other hand, if we retain only the terms of the first excited states in the 
summation, Eq.(36) gives the lower bound of Cq, namely 

~ 2 _ 

where Ri relates to the dipole transition amplitude z by 

R^ = 3\z\'^ with z = J d^ripl^^ir cos 9) i)o^o,o ■ (40) 

Since we know the value of Cq, Eqs.(38) and (39) will be used in turn to give the 
lower and the upper bounds of the strength of the underlying Coulombic force 
respectively, and which are 



f-C,iAE,) l<V<^^Ce(AEO (41) 

The numerical values of C and (3 of Eq.(29) for 7 = 1.5 are obtained by shrinking 
the energy range to be fitted. The results are 

7 = 1.5 , p = 0.0547 and C = 0.174 for [sc] (42) 
7 =1.5 , /? = 0.0541 and C = 0.172 for [ fc] (43) 

in the unit of the neutral pion mass, and the energy range of the fits are (0.6 ~ 34) 
MeV. and (0.6 ~ 39) MeV. for [sc] and [ fc ] respectively. Since the coefficient Ca 
of the asymptotic potential is 

C^ = C^ ^^^"~^^ with a = 27 + 3 , (44) 

Cq = 0.9933C" in our unit. If we use the mass difference of A(1232) and as the 
first excitation energy, namely AE = 2.18, then the lower bound of becomes 



> — ^ — • (45) 



r 



1 
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To go further we need the size of the composite particle, the radius of proton 
in particular. Although we do not have much information on the radius, we shall 
consider two cases. In the first case, is (0.7fm.)^, which comes from the nucleon 
form factor. In another case rf is set equal to /3 of Eqs.(42) and (43), because at 
the range r = ^/J3 the potential starts to deviate appreciablly from the asymptotic 
form of the van der Waals potential of V{R) ~ —Cq/R^ and therefore ^/J3 may be 
regarded as the size of the composite particle. The results are : 

> 3.27 ( from the form factor ) 

> 13.8 ( from ^ = /? ) . 

As it is expected, the underlying Coulombic force is super-strong, because the in- 
duced van der Waals force has already the order of magnitude of the strong inter- 
action. 

The coupling constants of the QED and the QCD are too small, they are 
(137.036)"^ and 0.32 respectively. On the other hand, the Coulombic force of the 
magnetic monopoles is an imortant candidate of the underlying Coulombic interac- 
tion. It is well-known from the charge quantization condition of Dirac that|jl9| 

* o 137.036 , , 

*e' = (46) 

Therefore of the magnetic charge is large enough to satisfies the condition of 
the lower bound, moreover it is not very large compared to the value of the lower 
bound. This situation is satisfactory, because the upper bound of Eq.(41) must have 
the value of the same order of magnitude as the lower bound. 

Since the Coulombic interaction between the magnetic monopoles is a suitable 
candedite of the underlying dynamics of the nuclear force, it is worthwhile to repeat 
briefly the magnetic monopole model of hadron0. Although it is overshadowed by 
the QCD, it still has virtue to be considered. The monopole model of hadron is 
essentially the quark model, in which the quarks bear the magnetic charges. Such 
a fundamental particle is often called dyon, since it is doublly charged, namely 
electrically and magnetically charged. Because of the superstrong Coulombic force 
between magnetic charges, the dyons form the bound states of magnetic charge zero, 
and such composite particles are identified with the hadrons. In a sense, hadrons 
are regarded as " magnetic atoms " in this model |0. Therefore it is not surprising to 
observe the strong van der Waals forces between hadrons. In this article we actually 
observed the strong van der Waals interaction in the proton-proton scattering. Since 
the van der Waals interaction is universal|2^, we may expect to observe such a 
force also in other scattering processes. However because of the lower precision 
of the data, the elimination of the two-pion exchange cut is necessary at least in 
the threshold region of the spectrum. After carrying out such eliminations, we can 
actually observe the long range force in the P-wave amplitudes of the ir-N and the 



TT-TT scatterings [21] 22 



20 



7 Remarks and Comments 



In this paper we searched for the long range force in the S-wave amphtude of the 
proton-proton scattering. In order to see the extra singularity of the long range force 
in holu) at u = as clearly as possible, we constructed a regular function KqIu) from 
the modified effective range function Xo(z/) of the proton-proton scattering, where 
the function KqIu) does not have the normal singularities in the neighborhood of 
1^ = 0. By separating the OPE part and then making the once subtracted function 
^once^j^^^ it becomes much easier for us to observe the extra singularity a.t u = 
when it exists. It is impressive to observe that K^^'^^ii') has a sharp cusp at = 0, 
and the form of which is close to the square root type, namely to Kq^'^'^{v) = 
(Cq + C"y/u) with positive C". The y/u singularity is exactly what is expected in 
the van der Waals interaction of the London type. Moreover the observed sign of 
C" indicates that the long range force is attractive. 

If we state it more precisely, when the spectral function has the form of (4m^, t) = 
nC't'^ exp{—pt), the chi-square fit to the multi-energy phase shift data[^ deter- 
mines the three parameters involved. They are 7 = 1.54, f] = 0.063 and C = 0.175 
in the unit of the neutral pion mass. It is interesting to see that the spectral function 
Af(4m^,t) has a peak at = 4.9 namely at 660MeV.. On the other hand, the lo- 
cation of the peak of the spectrum At{4:m'^, t)/t oi the once subtracted amplitude is 
y/t = 2.93 or 395MeV.. When we do not need the very precise description of the low 
energy phase shifts, the spectral function can be replaced by a 5-function located 
somewhere between 400 and 650MeV.. Such a 5- function may correspond to the 



fictitious (T-meson of the one-boson exchange model of the nuclear potential p3| p4 
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Since the power 7 of the threshold behavior of the spectral function At{Am?^ t) is 
close to that of the van der Waals force of the London type, we may assume that the 
ordinary mechanism, which produce the van der Waals force, is working. Then the 
strength of the underlying Coulombic force can be estimated. This is possible 
because there are the well-known inequalities on the coefficient Ce of the van der 
Waals potential |TB[, in particular the upper bound of which is obtained by using 
the strength and the radius ( )^/^ of the composite particle. Since we know 
the strength Cg of the van der Waals potential from the fitting, the inequality gives 
instead the lower bound of *e^, which depends on size of the composite state. For 
( y/'^ = 0.5, which is obtained by the measurement of the nucleon form factor, 
we obtain > 3.3, whereas for ( y/"^ = ^/J3 = 0.25, we obtain > 14. In any 
case, the underlying Coulombic force is super-strong. Among the interactions of the 
Coulombic form, the force between the magnetic monopoles satisfies this inequalty, 
because = 137.036/4 due to the charge quantization condition of Dirac[p!9| . 



Therefore the dyon model of hadron, in which the constituent particle dyon bears 
the magnetic charge as well as the electric charge, must be an important candidate 
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of the model of hadron[^. 

Finally, in order to confirm the long range force in the proton-proton scattering, 
it is desirable to observe directly the difference of the interference patterns by mea- 
suring precisely the angular distributions of the cross section in the low energy exper- 
iments. Our proposal is to observe the interference between the repulsive Coulomb 
force and the attractive long range interaction obtained in the present article. Since 
the interference pattern of the conventional short range foece and the Coulomb force 
is known, the difference can be predicted by using the values of parameters 7, 13 and 
C. It turns out that the most favorable energy is in Tiab = 25MeV. ~ AOMeV.. Since 
the relative deviation, namely A{d(j/dQ) ■ {da/dQ)~^, has a dip around 6c.m. = 10° 
with the depth ~ 0.3%, there exist two challenging requirements in the measure- 
ments. They are the accuracy and the domain of Oc.m.- However since the shape of 
the interference pattern is important, there is no severe conditions on the normal- 
ization of the cross section. Details will be published in a separate paper. 
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